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In this paper we investigate the speed of convergence of the 
fluctuations of a general class of Feynman-Kac particle approxima¬ 
tion models. We design an original approach based on new Berry- 
Esseen type estimates for abstract martingale sequences combined 
with original exponential concentration estimates of interacting pro¬ 
cesses. These results extend the corresponding statements in the clas¬ 
sical theory and apply to a class of branching and genealogical path- 
particle models arising in nonlinear filtering literature as well as in 
statistical physics and biology. 


1. Introduction. Feynman-Kac distribution flows and their particle in¬ 
terpretations arise in the modeling and the numerical solving of a vari¬ 
ety of problems including directed polymer simulations in biology and in¬ 
dustrial chemistry, nonlinear filtering in advanced signal processing and 
Bayesian statistics methodology, rare event estimation in telecommunica¬ 
tion and computer systems analysis as well as physics in the spectral analy¬ 
sis of Schrodinger operators and in the study of particle absorptions. Their 
asymptotic behavior as the size of the systems and/or the time parameter 
tend to infinity has been the subject of various research articles. For more 
details on both the theoretical and applied aspects of the topic we refer the 
reader to the review article [5] and references therein. 

To better connect this study with existing and related articles in the liter¬ 
ature we give a brief discussion on the fluctuation analysis of these models. 
The first “local” central limit theorems were presented in [1], These fluctu¬ 
ations were restricted to local sampling errors of an abstract class of genetic 
type particle model. This study was extended in [2] in the spirit of Shiga and 
Tanaka’s celebrated article [9] to particle and McKean path-measures. This 
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approach to fluctuations in path space was centered around Girsanov type 
change of measures techniques and a theorem of Dynkin and Mandelbaum on 
symmetric statistics [7]. This strategy entirely relies on appropriate regular¬ 
ity conditions on the Markov kernels which are not satisfied for genealogical 
tree evolution models as the ones described in [1]. Another drawback of this 
approach is that the description of resulting limiting variance is not explicit 
but expressed in term of the inverse of an L 2 integral operator. 

Donsker type theorems and an explicit computation of the limiting vari¬ 
ance in terms of Feynman-Kac semigroups were further developed in [4] in 
the context of particle density profile models. These explicit functional for¬ 
mulations were the starting point of a new approach to central limit theorems 
based on judicious martingale decompositions and Feynman-Kac semigroup 
techniques [3, 5]. 

The main objective of the current article is to complete and further extend 
these studies with the analysis of the speed of convergence of fluctuations. 

The article is organized as follows. In Section 1.1 we describe the Feynman- 
Kac and the particle models discussed in this article. In Section 1.2 we 
present our main results and specify the set of regularity conditions needed 
in the sequel. Section 2 is concerned with a precise Berry-Esseen type es¬ 
timate for abstract martingale sequences. In Section 3 we show how these 
martingale fluctuations apply to a sufficiently regular class of McKean par¬ 
ticle interpretations. 

We end this section with some rather standard and classical notation that 
will be of current use in the article. 

By M(E) we denote the set of all bounded and positive measures on a 
measurable space (E,£), by V(E) C JA(E) we denote the subset of prob¬ 
ability measures on ( E ,£) and by Bb{E) the Banach space of all bounded 
^-measurable functions f on E equipped with the uniform norm ||/|| = 
sup x eE |/(x)|. We also let Osc(E) C Bb(E) be the subset of all bounded 
measurable functions with oscillations osc(/) = sup( xy ) | f(x) — f(y)\ < 1. 

We finally recall that a bounded and positive integral operator Q from (E, £) 
into another measurable space (E',£') generates two operators, one acting 
on functions / £ Bb(E') and taking values in Bb(E), the other acting on 
measures n £ M.(E) into M.(E') and defined by 

Q(f)(x)= / Q(x,dx f )f(x f ), nQ{dx')= / fx(dx)Q(x,dx'). 

Je 1 Je 

To clarify the presentation we shall slightly abuse the notation, and we often 
write Q(f — Q{f)) 2 for the function x —> Q(f — Q(f)(x)) 2 (x). 

Finally, we shall use the letter c to denote any nonnegative and universal 
constant whose values may vary from line to line but which does not depend 
on the time parameter nor on the Feynman-Kac models. 
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1.1. Description of the models. We consider some collections of measur¬ 
able spaces (E n ,£ n ) ne n, of Markov transitions M n +i(x n , dx n +i) from E n to 
E n+ 1 , and bounded £ n -measurable and strictly positive functions G n on E n . 
We assume that the latter are chosen so that for any n G N we have 

(1) r n = sup (G n (x n )/G n (j/ n ))<oo. 

( x n ,Vn ) £ 

We associate to the pair ( G ni M n ) the Boltzmann-Gibbs transformation 
on V(E n ) and the mapping from V(E n ) into V(E n+ 1 ) given for any 

(x n -,hn) G (E n ,V{E n )) by 

(/In) {dXn) — G rl (x rl )/i rl ((ix n ,)//i n (G n ), 

4?n+l(/ln) = 'f , n(/ln)-M n -|-l. 

For any rj q G V(Eq) we denote by E, ?0 (-) the expectation operator with 
respect to the distribution of a Markov chain X n with initial distribu¬ 
tion 7/o and elementary transitions M n . We consider the distribution flow 
rj n G V(E n ), n G N, also called Feynman-Kac flow in the sequel, defined for 
any f n G Bb(E n ) by the Feynman-Kac formulae 


Vn(fn) =7n(/n)/7n(l) 


(2) 


With Tn (fn) = E 


' 7?0 


fn(X n ) [] G p (X p ) 

0 <p<n 


with the convention Yl 0 = 1. Using the multiplicative structure of the Feynman- 
Kac model and the Markov property, one readily checks that the flow ?/ n 
satisfies the nonlinear equation 


(3) Vn+1 — r/nKn+i^n , 

where (Kn + i i/ln ) ngN /ln6 p(g n ) is a nonunique collection of Markov transitions 
satisfying the compatibility condition 


(4) V7i G N V/r n G 7^(F/ n ) fi n K n+ i jPn — (/i n ). 

These collections of transitions are often called the McKean interpretations 
of (3). Notice that the compatibility relation (4) is satisfied if we take 

(5) E n +\ Prl ( X n , •) — £ n G n (x rl )M n _|_i ( X n , •) + (1 £ n G n (x n ))*& n - 1_1 (/In) 

for any nonnegative constant e n such that £ n G n (x n ) G [0,1]. We finally notice 
that the random variables X n may represent the path of an auxiliary Markov 
chain X' p from the origin up to time n and taking values in some Hausdorff 
topological spaces E' \ that is, we have 

(6) X n = (X' 0 ,...,X' n )eE n = (E' 0 x---xE' n ). 

For each iV>lwe denote by tiitv the mapping from the product space E N 
into V{E) which associates to each configuration x = ( x l )i <i<N G E N the 
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empirical measure rn/v(x) = J2i=i 5®*- The interacting particle system as¬ 
sociated to a given McKean interpretation is defined as a Markov chain 
taking values in the product spaces E^ with initial 
distribution ri^ N and elementary transitions 

N 

(7) Prob(^) G dx n |£>) = Ii K n,m N ^_\)^n-l\ dx n), 

i= 1 

where dx n = X i^^jydx^ stands for an infinitesimal neighborhood of the 
point x n = «)i <i< N G E%. 

Under appropriate regularity conditions on the McKean transitions ker¬ 
nels K n n n it is known that in some sense the particle measures 

Vn =m N (^ n N) ) 


converge as N tends to infinity to the desired distributions r] n . 

To illustrate this model we note that the particle interpretation of the 
Feynman-Kac flow associated to McKean transitions (5) forms a two-step 
selection/mutation genetic algorithm. The particular situation where e n = 0 
corresponds to a simple genetic model with an overlapping mutation/selec¬ 
tion transition. In the same vein the corresponding particle interpretation 
model of the Feynman-Kac path measures associated to the chain (6) forms 
is a genetic type algorithm taking values in path space. Note that in this 
situation the path-particles have the form 


S n 


( JN,i) JN,x) 

v>0,n JSl } 7i 



(E' 0 x ■ • ■ x E' n ). 


In addition, if the potential functions only depend on terminal values in the 
sense that G n (x' 0 ,... ,x' n ) = G' n (x ' n ) for some potential function G' n on E' n , 
then the resulting path-particle model can be interpreted as a genealogical 
tree evolution model. 

As traditionally, to clarify the presentation we slightly abuse the nota¬ 
tion, by suppressing the size index N, and we write (m(a;),£ n ,£^) instead of 
(m N (x),& N) 


1.2. Statement of some results. For any sequence of JF^-adapted random 
variables Z^ defined on some filtered probability spaces (£l N , (F^) n >o,F N ), 
we denote by Athe difference process A = Zff — Z^ f _ 1 , with the con¬ 
vention AZq = Zq for n = 0. If AM,f is a given tF *-martingale difference, 
then we denote by the jF^-martingale defined by = J2p =o AM^. 
We recall that its increasing process ( M N ) n is given by 

n 

(M N )n = J2 k n [(AM») 2 \F?_ 1 }, 

p=0 
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with the convention FYi = {0, for p = 0. It is also convenient to in¬ 
troduce the increasing process C^ = N(M N ) n of the normalized martingale 
Ln = VNM£. 

The example we have in mind is the situation where o is the 

natural filtration associated to the particle model (7) and the JP v -rriartingale 
difference A (f n ), with f n €B b (E n ), is given by the particle nth 
sampling error 

(8) A MZ(f n ) = rfi(f n )-rfi_ 1 K nir gJf n ) 

with the convention 7]Y\K 0t) n^ = r/o for n = 0. In this situation the increas¬ 
ing processes (A i N (f)) n and (/) = (£ N (f))n of the corresponding mar¬ 
tingales Mn (/) and C% (/) = \/N (/) are connected by the formula 


(9) C%(f) = N(M N (f)) n = K p ^Jf P 


-K 


P’% 


p =o 


Jf P )Y 


Our first main result concerns a Berry-Esseen theorem for an abstract class 
of martingale sequences under the following set of conditions: 

(HI) For any n > 0 there exist some constants oi(n) < oo and 0 < ci(n) < 1 
such that for any n > 0 and A 3 < ci(n)A 1 / 2 we have, almost surely, 

]E[e iXN^AM»+(*mAC» l:F N _^ _ 1| < ai(n)A 3 /iV l/2_ 


(H2) For any n > 0 there exists some finite constant 02 (n) < oo such that 
for any N > 1, A > 0 and n > 0, 

| E [ e UA 1/2 M£ r ]| < E j e -(A 2 /2)AC^] e A 3 a 2 (rx)/IV 1 /2^ 


(H3) There exists a nonnegative and strictly increasing deterministic pro¬ 
cess C = (C' n ) n >o as well as some finite constants 0 < a$(n) < 00 such 
that for any e > 0 we have 

e ^ivV2|ac^-ac„|] < (1 +eo 3 (ra))e e2a 3( n ). 


Theorem 1.1. Let M N = (M^) n > 0 be a sequence of F N -martingales 
satisfying conditions (H1)-(H3) for some nonnegative and strictly increas¬ 
ing process C n . We let F^ , respectively, F n , be the distribution function of 
the random variable , respectively, the one of a centered Gaussian 

random variable with variance C n . Then for any n > 0 we have 

limsup7V 1/2 ||F ? f — F n \\ < 00 . 
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Theorem 1.1 does not apply directly to the particle martingale sequence 
introduced in (8). The first two conditions (HI) and (H2) are rather stan¬ 
dard. They can be checked for any kind of any McKean interpretation model 
using simple and rather standard asymptotic expansions of characteristic 
functions. The third condition is an exponential continuity condition of the 
increasing processes introduced in (9). Next we provide a sufficient regular¬ 
ity condition which can be easily checked in various McKean interpretation 
models. If we set for any fi n E V(E n ) 

/■^(O) = {he Osc {E n ): n n {h) = 0}, 
then this condition reads 

(H) There exists a collection of uniformly bounded positive measures T n j 
and T( l+1 j on the sets ^ 1 (0) and r/^, (0) and indexed by n E N and 
/ E Osc(£' n+ i) and such that 

\\K n+ 1 ^(f)-K n+ l,r, n (f)\\ 

< j \fin(h)\T nJ (dh) + J\$ n+1 (n n )(h)\r' nJ (dh). 

When condition (H) is met we denote by T the supremum of the total mass 
quantities V n ^(1) and T n j(l). 

Note that (H) is related to some Lipschitz type regularity of the increasing 
process and it is clearly met for the McKean transitions given in (5), since 
we have in this case 


Kfi+l ,/inif) ^n+l,T)r,(/) — (1 £nG n )[*&n+l(Hn) ^n+l {Vn)] (/) ■ 

Thus, in this situation, we have that (H) is met with T n j = 0 and T^ j- = 
§h with h = [f — rj n+ i(f)] so that T n j( 1) = 0 and T^y(l) = 1. When the 
parameter e n = e n (/r n ) in (5) depends on the index measure /r n , we also find 
that (H) is met as soon as we have 

\£n(fJ’n) E-nijln) \ ^ J \l^n (^) l-A-n (dfl) 
for some uniformly bounded positive measures A n on ?7“ 1 (0). 

Remark 1.2. The above considerations show that condition (H3) is in 
fact easily verified in most of the classical applications of McKean models, 
and in particular in the case of nonlinear filtering, for which we refer to [5], 
for sake of conciseness. 

To check the exponential estimates stated in condition (H3) we shall use a 
refined version of Burkholder’s inequality recently presented by the first au¬ 
thor with Miclo and Viens in [6]. Roughly speaking, these sharp L p -estimates 
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combined with some judicious error decompositions lead to the desired ex¬ 
ponential concentrations estimates for the normalized sampling error mar¬ 
tingale (/) defined by (9) with the limiting increasing process 

n 

( 10 ) Cn{f) = '^Vp-lKp, r lp-l(fp — Kp, r i p - 1 (fp)) ■ 

p=0 

Observe that, even if this strategy led to the desired Berry-Esseen esti¬ 
mates on (/), we would still need to transfer these rates of fluctuations 
to the random field sequence defined by — rj n ). One of the most 

elegant approaches is probably to follow the semigroup techniques and the 
martingale decompositions developed in [5]. To describe these decomposi¬ 
tions with some precision we let Q V)U be the linear Feynman-Kac semigroup 
associated to the flow 7 n . To be more specific, we define the semigroup Q P)U 
by the relation 


7 n — 7 pQ 


p,n, 


and we associate to Q p n a “normalized” semigroup Q p n , defined for f n £ Bb{E n ) 
by 


(ll) 


Qp,n(fn ) — 


Qp,n{fn ) _ 7p(l) 

VpQp,n{ 1) 7n(l) 


Qp,n{fn) • 


If we let {W pn {f n )) p < n , f n £ t3b(E n ), be the random field sequence defined 
by 


(12) W p N n (f n ) = VN(r,» -r, p )(f p , n ) 


with f P) n = Q Pj n(fn ~ Vnfn), 


then we have the Doob type decomposition 

K n {fn) = KnUn) + C^ n {f n ), 

with the predictable and martingale sequences given by 

A<„(/„) = Vn[ 1 - 

X [^p(Vp—l){fp,n) — ®p{Vp-l)(fp,n)], 

A £% n (fn) = VN[Vp(fp,n ) - 

The above decomposition is now more or less standard. For the convenience 
of the reader its proof is housed in the Appendix. 

Intuitively speaking, we see from the quadratic structure of the predictable 
term that it should not influence the fluctuation rate. We will make precise 
this observation with a Stein type approximation lemma and we will prove 
the following theorem: 
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Theorem 1.3. Let f n E Bb{E n ) and let W^ n be the quantity defined 
by (12). For any McKean interpretation model satisfying condition (H), we 
have 



for any f n E Bb(E n ) and n > 0 with 


n 



2. An estimate for martingale sequences. The central limit theorem for 
sequences of random variables is usually obtained by convergence of char¬ 
acteristic functions. Unsurprisingly, the natural question of determining the 
speed of convergence in the CLT can also be handled through characteristic 
functions considerations. The formalization of this idea is due to Berry and 
Esseen, and can be summarized in the following theorem: 

Theorem 2.1 (Berry-Esseen). Let (Fi. F 2 ) be a pair of distribution 
functions with characteristic functions (/i,/ 2 )- Also assume that F 2 has a 
derivative with ||^.|| < 00. Then for any a> 0 we have 



In this section, we will try to apply this theorem to a sequence of mar¬ 
tingales satisfying the general set of hypotheses (H1)-(H3) in order to get a 
sharp asymptotic result for its convergence towards a Gaussian martingale. 
In order to prepare for the proof of Theorem 1.1 we start with the following 
key technical lemma. 

Lemma 2.2. Suppose we are given a sequence of martingales = 
0 with respect to some filtrations , satisfying the conditions (Hl)- 
(H3). Then, for any n> 0, there exist a finite constant a(n ) < 00 , a posi¬ 
tive constant b(n ) and some N(n ) > 1 such that for any N > N(n) and 
0 < A < b(n)N 1 / 2 , 


| E[e *A _ e -(A»/2)C B| < a(n)e -(A 2 /4)AC„^(l + A) ^ 


Since the proof of Theorem 1.1 is a simple consequence of the above lemma 
we have chosen to give it first. 
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Proof of Theorem 1.1. By Theorem 2.1 and Lemma 2.2 we have, 
for any N > N (re), 


iV 1 / 2 ||F n Ar -F n || 

< ^ [ b{n)Nl/2 e~^^ AC -X(l + A) dX 
vr Jo 


+ 


24 


b(n)(2eTr 3 C n ) 1 / 2 


< 


2 a(n) 


IT 


poo 

[ e -(A 2 /4 )AC„ A (i + A ) dA + 

Jo 


24 


b(n)d /2 ' 


for some N(n) > 1 and some finite positive constant 0 < b(n) < oo. Invoking 
the fact that A C n >0, this ends the proof of the theorem. □ 

We now come to the proof of Lemma 2.2. 

Proof of Lemma 2.2. Let 1 ^ be the function defined for any A > 0 
by 

/^(A)=E[e iAArl/2M ^ 2 /2)C' n] _ 1) 

and notice that 

(13) E[e iXNl/2M "} - e -( A2 /2)C n = e -(A 2 / 2 )C n/ A(A)_ 

Furthermore, we have the easily verified recursive equations 

InW-In-lW 

= E [ e iAI vl/2 M~_ i +(A 2 / 2) C n _ i ( E [ e aA 1/2 AM^+(A 2 / 2 )AC„| F jV j 
and hence 

I»(X)-I»_ l (X)=A(X) + B(X), 

with 

A(X) =E [ e iAiV 1/2 M^_ 1 +(A 2 / 2 )C n - 1 ( e (A 2 / 2 )(AC„-AC^) _ ^ 

and 

B(X) =E[e iXNl/2M ^-i + ^ 2/2)Cn - 1 (e^ 2/2){ACn ~ AC ^ ) ) 
x (E[6 iAJVl/2AM »+( Aa / 2 > AO » |^_i] - 1)]. 

Using this we obtain 

\C(X)-I^_ 1 (X)\<e^^ c ^(A 1 (X) + B 1 (X)), 

where 

A 1 (A)=E[e( A2 / 2 ^ AC "- AC »l-l] 
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and 


£?!(A) =E[|E[e aJvl/2AM " Ar+ ( A2 / 2 ) A ^|^_ 1 ] - i| e (A 2 / 2 )(AC n -AC^) ] _ 
Now, under conditions (HI) and (H3) applied for e = A /2 > we fi R d that 
|J*(A)-J*i(A)| 


< e 


(A 2 /2)C n _i 


ai {n)\ 3 f i | A 2 a 3 (n) \ (A 4 /4jV)a 2 (n) 


IV 1 / 2 


21V 1 / 2 


= e (A 2 /2)C n _i 


/ A 2 a 3 (w) A (A 4 /4JV)a 2 (n) _ 

+ V 21V 1 / 2 y 

Q l( n ) A ‘ 3 ( 1 + A2 «3(n) \ p(A 4/4JV)a 2 (n) 


IV 1 / 2 


21V 1 / 2 


. r „(A 4 /4N)a 2 (n) _ ^ (A 4 /4A)a 2 (n)' 

+ l 21V 1 / 2 

for any 0 < A 3 < ci(n)lV 1 / 2 . Since for these pairs of parameters (A, IV) we 
have A 2 < IV 1 / 2 (and therefore A 4 < IV), we find that 

1V 1/2 |^(A) - lA_i(A)I < d(n) e ( A2 / 2 ) c "- 1 A 2 (l + A), 

for some finite constant d(n), whose value only depends on aj(n), i = 1,3, 
and such that 


d(n) < ce a2( - n ^ 4 (l V ai(n) V a 3 (n))^ 


If we set 


c*(ra) = A c i(?>) (- !) and d *(n) = V d (p)> 

p=0 p=0 

then for any 0 < p < n and any 0 < A 3 < c^r^lV 1 / 2 , we have that 

lV 1 / 2 |lA(A) -AA( A )I <^We (A 2 / 2 )C "- 1 A 2 ( 1 + A). 

It is now easily verified from these estimates that 

!V 1 / 2 |lA(A)| < (n + l)cf(n)e( A 2 / 2 ) c - 1 A 2 (l + A), 

from which we conclude that, for any 0 < A 3 < c*(n)lV 1 / 2 , 

(14) \E[e lXNl/2M "} - e~^ x2/2)Cn \ < ( n +l)d*(n)e~^ 

IV 1 / 2 

On the other hand, we have, for any pair (A, IV), 

(15) | E [ e ^ 1/2 M^] _ e -^/2)C n \ < |E[e*AN 1 / 2 M^j| _j_ g- 


(A 2 /2 )C n 
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and under condition (H2), 

j | < E j e -(A 2 /2)A C» ] e A 3 a 2 (n)A-V2 _ 

Using again (H3) we also find that 


|E[e' 


iAA 1 / 2 M^i 


< -(A 2 /2)Ac„ ( 1 , A 2 q 3 (n) \ ( A 4 a 2 (ri)) /( 4jV ) (A3 ffl2 ( n) ) /jV i/ 2 

- ^ + 21V1/2 ; e 

A 2 a 3 (n)^ A 2 / 


= 1 + 


21V 1 / 2 


exp 


2 V 


ACL — 


N V 2 


2a 2 (n) + a 3 (n) 


A 


21VV2 


Recall that A C n > 0, and observe that for any pair (A, IV) such that 
A < c*(n)N 1 / 2 with c*(n) = [2ag 2 (n) A (2 _1 AC n (l + 2a 2 (n)) -1 )], 


we have 


jyd /2 ( 2a 2(") + ^(2a 2 (n) + 1) 

< ACL 
“ 2 


This yields 


(16) 


|E[e : 




< 1V 


a 3 (n) 


i + A2 ^ p -A 2 AC' n /4 

i+ ivv2; e 


and hence, by (15), and for any A < c*(n)lV 1//2 , we find 


(17) |E[e lAJVl/2M "] - e -( A2 / 2 )C„ , < e -(A 2 /4)AC„ (2 v 03 ( n)) (i + . 


To take the final step we observe that for any 

A > c*(n)/c*(n) 3 and cy 3 (n)A 1 / 6 < A < c*(n)A 1//2 , 
we have 1 = c*(n)/c*(n) < c“ 1 (n)A 3 /lV 1 / 2 , and by (17), 

(18) \E[e iXNl/2M "} - e - (A2/2)c "| <c; 1 (n)(2Va 3 (n))^y^e- A2AC " /4 . 

In conjunction with (14) we conclude that for any N > N(n) = c*(n)/c*(n) 3 
and any A < c*(n)N 1 / 2 , 


| E[e IAiW 2 M^ ] _ g — (A 2 /4)C n | < a(n )_^(l + A)e -(A 2 /4)AC n) 

TV 1 / 2 

with a(n) = [(n + l)d*(ra)] V [c“ 1 (n)(2 V a 3 (n))]. This ends the proof of the 
lemma. □ 
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3. Application to interacting processes. In this section, we prove that 
Theorem 1.1 can be applied to our particle approximations. We shall go 
through a series of preliminary results leading to the proof of Theorem 1.3. 

The first step is of course to provide some exponential estimates for the 
particle density profiles. In the next pivotal lemma we describe an original 
exponential concentration result in terms of the following pair of parameters: 

P(Pp,n) = Slip ||/p jrl (Xp, •) Pp,n{ypt •) 11tv, 

(xp,Vp)&El 

(19) 

fp,n — SUp Qp,n(l) i%p)/Qp,n (1) (?/p); 

where \\p — v || tv = sup^ g £ |/r(A) — v{A)\ represents the total variation dis¬ 
tance between probabilities and Pp )U denotes the Markov transition from E p 
into E n defined by 

Pp,n(%p > dx n ) = Qp^niXpi dXfi)/Qp : n(%pi E n ). 


Lemma 3.1. For any McKean model we have for every n > 0, f n G 
Osci (E n ) and e > 0, 

E [ e ^ 1/2 K(/n)-»?n(/n)|] < (i + e 2“ 1/2 6(n))e (e6(ri))2/2 
for some finite constant b(n) such that h{n) < 2J2q = o r q,n/3(Pq, n ) ■ 

Remark 3.2. The quantities (r P:Tl , /3(P Pjn )) play an important role in 
the asymptotic and long time behavior of Feynman-Kac particle approxi¬ 
mation models. The above lemma combined with the semigroups approach 
developed in [5] readily yields uniform exponential concentration properties. 
To be more specific, let us suppose that r = \/ n r n < oo. Also assume that 
there exist some integer parameter m > 1 and some p G (0, 1 ] such that for 
any (x , y) G E^, A G £ n+m and n > 0, 

M n ^ 4 ) pM n ,n-{-m(yi Af 

where M n ^ n+m = (M n+ 1 ... M n+m ) stands for the composition of the Markov 
kernels M p from p = (n + 1) to p = (n + m). In this situation, following the 
arguments given in [5] one proves that 

r n ,n +m <r m /p and fi(P n , n+m ) < (1 - . 

Furthermore, the constants b(n) in Lemma 3.1 can be chosen such that 
\J n b(n ) < 2mr 2m ~ 1 / / o 3 . 
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The proof of Lemma 3.1 being rather technical, it is housed in the Ap¬ 
pendix (see Lemma A. 3). One consequence of Lemma 3.1 is the following 
central estimate. 

Lemma 3.3. Suppose the McKean interpretation model satisfies con¬ 
dition (H) for some finite constant T < oo. In this situation, the martin¬ 
gale (/) defined by (8) satisfies conditions (H1)-(H3) for some universal 
constants 

(ai(n),a 2 (n)) = (ai,a 2 ) 

with the nonnegative increasing process C n (f) defined at (10), as soon as 
n—>C n (f) is strictly increasing. In addition, the constant a^in) in (H3) 
can be chosen such that, for any n > 0, 

<? 

0 < 03 ( 71 ) < 4a/ 2(1 + T) sup 

q=n,n-l p=0 

Furthermore, when the regularity conditions stated in Remark 3.2 are 
met for some triplet ( m,r,p), the constant a^^n) can be chosen such that 
0 < V n a 3( n ) — 8v / 2mr 2m ^ 1 (l + T)/p 3 . 

The second step will be to get rid of the predictable term defined by (12), 
with the help of the following lemma: 

Lemma 3.4. Let Fz be the distribution function associated to a real¬ 
valued random variable Z, and let W be a centered Gaussian random variable 
with unit variance. For any pair of random variables (X,Y) we have 

(20) || F x +y -F w \\< || F X - F W \\+4E(\XY\) + 4E(|y|). 

Lemma 3.4 can be proved using the Stein approach to fluctuations and 
it can be found, for instance, as Lemma 1.3, Chapter 11 in [8]. Since the 
proof of Theorem 1.3 is now a simple consequence of Lemmas 3.3 and 3.4, 
we postpone the proof of Lemma 3.3 and start with the 

Proof of Theorem 1.3. Throughout this proof, b(n ) will stand for 
a positive constant that can change from line to line. We first notice that 

W% n = VN(r,% - Vn ), 

and by (12) we have W pn {f n ) = B pn ff n ) + C pn {f n ). Let us show now that 
the main term in the fluctuations of the c.d.f. of W pn (f n ) is due to C p n (f n ). 
Indeed, it is easily checked that 

E 1/2 (|<n(/n)I") V nK n (fn)\) < 


(21) 
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By definition of the martingale term, it is also easily seen that 

(22) JE 1/2 (I C^ n {f n )\ 2 )<b{n)(C n (f)) 1/2 . 

Set now 

X = £n,n(fn)/(C n (f)) 1/2 and Y = < n (/ n ) / (C n (/)) 1/2 . 

The estimates (21) and (22) yield 

EG^I) < ^E 1 / 2 [|< n (/ n )| 2 ]E 1 / 2 [|^ n (/ n )| 2 ] < 

and E(|T|) < b(n)/N 1 / 2 . Hence, applying Lemma 3.4, the proof of Theo¬ 
rem 1.3 can be reduced to show that 


sup 

WElL 


HC-n (/»,») < «(C„(/)) 1/2 ) - e^' 2 dv 


< 


b{n) 
N 1 / 2 ' 


This last estimate is now a direct consequence of Lemma 3.3 and Theo¬ 
rem 1.1. 

□ 


We now come to: 

Proof of Lemma 3.3. Let us first check that the regularity condi¬ 
tion (H3) is satisfied. Since we have 

ACn(f) = Vn-1 [Kn^Atfn - K n ^_J n ) 2 )} 

= $ n (?7 n -i)(/ 2 ) - rj n -i{{K ntVn _J n ) 2 ), 

we easily prove that 

IA C(/) - AC n (f)\ < 2(|<h n (r^_ 1 )(M + l^-i^-i)! 

+ \\K n r) N (f n ) — K n Vn _ 1 (f n )\\), 

’ 'n— 1 

with the pair of functions (h n ,h , n _ 1 ) £ (Osc (E n ) x Osc (E n _i)) defined by 

K = (fn ~ %(/n))/2, 

h n -l = ((Kn^n-ifn) ~ r ln-l((Kn, r i n _ 1 fn) ))/ 2 - 
On the other hand, under condition (H), we have that 

(fn) ~ K n ,r)„__ i(/n)|| 

< J\Vn-l(h)\T n -lJ n (dh) + J\*n(Vn-l)(h)K_ 1Jn (dh\ 
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from which we find that 

\AC%(f)-AC n (f)\ 

<2 (/ I^Wlfn-i^+Ed 

with 

f n -i = r„_ij n + 5^^ and f^i = + Sh n - 

Applying Jensen’s inequality, we get that for any e > 0, 
E [ e £A 1 /2|AC^(/)-AC n (/)|] 

< E[ e 2e7Vl/2{ /l^-iW|r»-i (dh)+f |^(fc)|r;_ 1 (dfc)}j_ 

Now, applying the Cauchy-Schwarz inequality, we obtain 
E [ e eA 1 /2|AC^(/)-AC„(/)|] 

< E 1 / 2 [e 4eArl/2 ./' ieiWI?-iW] E i/2 [e 4 £ ArVV |??^(M|r^ l _ 1 (d/i)j 
If we set r = r + 1, then using again Jensen’s inequality, we find 
E reA 1 /2|AC(/)-AC rl (/)h2 < f E r4 £ AV 2 f |r ? ^_ 1 (Mil ^-1 (dh) 

[ J ~J J r„-i(i) 

x f E[ c 4£jvV2?l ^' 

J K-iM ’ 

from which we get 

E[e£ AV2|AC(/)-AC„(/)l ] < sup E(e 4sAV2f| <W | ) _ 

h(E:Osc(Ep) ,p=n,n —1 

Using Lemma A.3, we conclude that 

E [e £ivl /2| A C(/)-AC n (/)|j < + £a3 ( n ))e e2a iW 

for some finite constant 03 (n) such that 

_ <? 

03 (n) < 4\/2T sup 

g=n,n—1 p_ 0 

To prove that (H2) is met, we first recall that 

| E [ e iANV2.M^ (/) ]| <E[|E[e aJV i /2 A - M « ( / ) |F^ 1 ]|]. 


(23) 
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Then we use a standard symmetrization technique: given the particle model 
up to time p < n — 1, we let fj* be an auxiliary independent copy of . 
In other words, fj is the empirical measure associated to an independent 
copy of the configuration of the system at time n. With some obvious 
abuse of notation, we readily check that 


| E [gUA 1 / 2 | pN_ j |2 = E [giAJV 1 / 2 [A (J)-AM %(/)] | pN_ 


n— lJ) 


where AM n (/) = [ 77 " (f n ) - $n(Vn-i)(fn)]- We deduce from this that 


|E[e 


i\N 1 / 2 AM%U)\F N u2 


N 

ei]i 2 =n^ 

i=i 






E[ e *(A/A rl ' /2 )[/ n (^)-/ n (^)]|j E Ar j =E 


'j^N 
S’ n—1 


Since the random variables [/„,(££) ~ fn(Cl)] and -[/«(££) - fn(£h)] have the 
same law, their characteristic functions are real, and we have 

cos (^72 [/**(&) - fniti)]' 

Using now the elementary inequalities 

cosu < 1 — u 2 /2 + |u| 3 /3!, l + u<e“, \u — u| 3 < 4(|u| 3 + |u| 3 ), 
we get that 

]E[g*(A/-N 1 / 2 )[/„(^)-/„(^)]|^ r A j 

\2 „\3 

< 1 - ]y K n,^_Sfn - K^N^fn)) 2 (&) + 

< e-^^K^N ^fn - K^N ^fn)) 2 ^) + 

Multiplying over j, we obtain 

m e i^ 1/2 ^ M n(f)\^_^2 < g—A 2 AC^ (f)+c\ 3 /N 1 / 2 ^ 

and by (23) we conclude that condition (H2) is met with 02 ( 71 ) = c/2. 

We now come to the proof of (HI). By definition of the particle model 
associated to a given collection of transitions K n ^ we have 

g^AJV 1 / 2 AM% (/)+(A 2 /2)AC*(/) | jrN_ j 
N 

= n[^e 1 ( el(A/ ^ 1/2)/ ' +(A2/(2 ^ ))AC " (/) )K^-i)> 

3 =1 

with the random function // = (/„ — If jv (/n))(6!_i)- Using the elemen- 

’ ''n — 1 

tary inequality 


e 2 - 1 + z + 


< e l 2 l. 


3! ’ 
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it is easily seen that, for any A < IV 1 / 2 , we have 
e i(A/IV 1 /2)/^ + ( A 2/(27V))AC (/) 

= 1 + (/) - Th 2 ]+<1 </), 

with \rJJ 1 (f)\ < cA 3 IV _3 / 2 . This clearly implies that, for any A < IV 1 / 2 , 
[*-„,» _ (e i(V-V'.,/> + (AV(2 I V))AC»(,) )]K J ) 


= 1 + ^[AC(/) - ^_ 1 (/^) 2 (^-i)] + < 2 (/), 

with |r , ^ 2 (/)| A cA 3 IV~ 3 / 2 . It is now convenient to notice that for any A < 
IV 1 / 2 , 


A[AC"(/) - K, 


N \ 
n ’Vn-l ' 


(/h 2 «Ti)]+<2(/) 


< 


cA 

JV 1 / 2 ' 


On the other hand, for any |z| < 1/2 and with the principal value of the 
logarithm we recall that 


/ ^ U I ^ L 

log (1 + z) = z — / -- du = z — z 2 / -- dt. 

J o 1 + u Jo 1 + tz 

Since for any \z\ <1/2 and t £ [0,1] we have |1 + > 1/2, we find that 

for any \z\ < 1/2 we have | log (1 + z) — z\ < \z 2 \. The previous computations 
show that there exists some universal constant Co £ (0,1) such that for any 
A < coiV 1 / 2 we have 


1 og^,<.,(= i(VNl/2,/ " +(XI/<2N)|iC - w )«hi: 


= A[a C(/) - (/j) 2 ^,)] + < 3 (/), 


with |r^ 3 (/)| < cA 3 /IV 3 / 2 . Summing over j, we see that for any A < co-ZV 1 / 2 , 

< cA 3 /IV 1/2 . 


A 


SlogK„, ei ( e «V'V' 2 )/i+(» a /(21V»A C »( / , )K j_ i; 

i=i 

Finally, using the elementary inequality \e z — 1| < |z|e^, we conclude that, 
for any A < Co-ZV 1 / 2 , 

| E [ e iAA 1 / 2 AA(^(/)+(A 2 /2)AC(/)| JZ rA j _ 1 | < c A ' j ^cX 3 /N 1 / 2 . 

N 1 / 2 

This readily implies that there exists some universal positive constant c\ 
such that, for any A 3 < c\N 3 / 2 , we have 

| E[e iAAV2AAf^(/) + (AV2)AC(/)|_ F A ^ < c _^_ 
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which proves that condition (HI) is met with a\(ri) = c and ci(n) = 1. □ 

APPENDIX 

A.l. Doob type decompositions. 


Proposition A.l ([5]). Let ( Q p , n )p<n be the semigroup defined at (11). 
For f n G Bb(E n ) and p<n we set f p ^ n = Q p ^ n {f n — Vnfn)- Then we have the 
following decomposition: 

Vp(f P ,n)=A^ n (fn) + M^ n (fn), 

with the predictable and martingale sequences A pn (f n ) and Ai^ n (f n ) given 
by 

(A.l) <„(/„) = til - <- 1 (Q,- 1 ,(<-,)(/,.«), 

q=l 

(A.2) (4«) - 

q=0 

with the usual convention <I>o( r ?-i) = Vo- 


Proof. Note that for any (p n G B^En) we have the decomposition 

_ _ p 

Vp {Qp,nTn) — Vo (Qo^Tn) = 7 ^ $qi 

q =1 

with 5 q = Vq (Qq : nTn) ~ Vq-1 (Qq-l^Vn) • Choose llOW ip n = fn~ Vnfn■ For 
q<p, we have, by definition of fg tn , 

$q = Vq ( fq,n ) “ Vq-l(Q q-l,nVn) = U 1 +U 2 , 

with 


Ul=V^(fq,n)-[%(V^lWq,n), 

U2 = [$q(Vq-l)\(fq,n) ~ Vq-l(Qq-l,nVn)- 
In order to show (12), it is thus enough to verify that 

( A -3) Vq-l(Q q -l,nTn) = Vq-l(Qq-l,ql)[$q(Vq-l)](fq,n)- 

However, we have 


7a .. _7?-i(l) r> .. _7 q-i(1)^^ .. A 

^iq—l.n^Pn — 7 (1) — 7 (1) 

_ TViUln (q 09 ) — 7 g-i(l) q / £ \ 

— ^ ^<lV°tq,nFn) — ^ ^q{Jq,n), 
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and hence 


(A.4) 


Vq-liQq-l^n) = ^ ^ r]q-l{Qq{fq,n)) 

Iq-lWVq-liGq-l) , 

7,(1) 1 9 




On the other hand, for q > 1 and x g _i € E q _ i, 



and 



which yields 



Plugging this last equality into (A.4), we get (A.3), and hence (12). The 
martingale property oi Mp n (f n ) is readily checked. □ 

A.2. Some asymptotic estimates. The next lemma provides a refined 
version of Burkholder type inequalities for independent sequences of random 
variables. 

Lemma A. 2 ([6]). Let m(X) = be the N-empirical measure 

associated to a collection of independent random variables X 1 , with respec¬ 
tive distributions /r* on some measurable space (E,£). For any sequence of 
£-measurable functions hi such that Hi(hi) = 0 and & 2 (h) = JffLi osc 2 {hi) < 
oo, we have for any integer p> 1 


y/NE(\m(X)(h)\ p ) 1 ^ p < d(p) 1/,p cr(h) 


(A.5) 


with the sequence of finite constants (d(ri )) n >o defined for any n > 1 by the 
formulae 


(A.6) d{2ri) = (2n) n 2 n and d.(9,n. — 1) = 1 / 2 \ 


y/n - 1/2 


The extension of Lemma A.2 to the interacting particle measures r}^ and 
the Feynman-Kac flow rj n is the following 
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Lemma A. 3. Let (d(p)) p > ± be the sequence introduced in (A. 6). For any 
McKean interpretation model and for any n > 0, p > 1, f n £ Osc (E n ) and 
e > 0, we have 

E(|[< - Vn\(fn) \ P ) 1/P < d(p)^b(n)/VN, 

E[ e e* 1/a \rfi(Jn)-TlnUn)\] < (1 +e2 

for some finite constant b(n ) such that b(n ) < ^J2q=o r q,nP(Pq,n)- 


Proof. The proof is based on the following decomposition: 
(A.7) 


Vn - Vn = '52[&q, n ( r ff) ~ ^qA^qiVq- 1 ))]- 
q =o 


We introduce the random potential functions 


G% n :x q €E q ^GZ n {x q ) = 


G 


q,n 


W-l )(G q ,n) 


£ (0, oo) 


and the random bounded operators P qn from Bb(E n ) into Bb{E q ) defined 
for any (f n ,x q ) £ ( B b (E n ) x E q ) by 

Pq!n(fn)( x q) = Pq,n(fn ~ ^ q,n(^ q(jlq-l )) (fn)) { x q) 

= f ( Pq,nf{ x q ) “ P^nf^q^G^yq^q^^dyq). 

We associate to the pair {G qn ,P qn ) the random bounded and integral oper¬ 
ator Q qn from Bb{E n ) into B b {E q ) defined for any (f n ,x q ) £ ( Bb(E n ) x E q ) 
by 


Qq,n(fn)( X q) = G q,n( X 1 ) X P q,n (fn) ( x q) ■ 

Each “local” term in (A.7) can be expressed in terms of Q q [ n as follows. For 
any q < n and f n £ Bb(E n ) with osc (f n ) <1 we have 


<1>,,»«)([/» - {<-!))(/»)]) 


N 


iq j\u n qa i \ q\ iq 

1 




^(Gq,nPq,n[fn ~ ))(/,)]) 


= 11qQq,n(fn) 

VqQq,n (!) ’ 

By construction we also observe that 

^q(Vq-l){Gq )n ) = 1 and <S>q(Vq-l)(Qq,n(fn)) = 0- 
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The above considerations easily yield the decomposition 

1 


$q,n(Vq) ~ ®q,n(®q(Vq- 1 )) = M <rN S ^q ~ $ q(Vq-l)]Qq,n- 


<(G» n ) 

Using the properties of the Dobrushin contraction coefficient, we also have 

\\P£ n (fn)\\<Oac(P q , n f)<P(Pq, n ), 

\\Qq,n(fn)\\ < \\G% n || \\P£ n (/„) || < \\G^ n \\P(P q , n ), 
and from these estimates, we readily prove the inequality 


PwLf) - ))](/„)! 


N 


< r q , n P(P q , n )\[rff ~ <S>q{r)q-l)}Qq,n(fn)\, 


N 


tN 


-ttN 


with Q q ] n (fn ) = Qq,n(fn)/\\Qq,n(fn)\\- Now, using Lemma A.2, we check that 
for any p > 1 we have 

VNE(\[rfi - %(V^l)]Qq,n(fn)\ P K-l) 1/P < 2 d( P ) 1/p , 

with the sequence of finite constants d(p) introduced in (A. 6). This ends the 
proof of the first assertion. The L n -inequalities stated in Lemma A.2 clearly 
implies that, for any e > 0, 


E(e‘ 


?\Vn _ 


= E 

n>0 


~2n 


(2 n)\ 


E (\Vn (fn)-Vn(fn) I 


2 n\ 


r 2n+l 


+ n? 0 ^ + 1 ) ! 


n\^(fu)-Vn(fn)\ 2n+1 ) 


<E 

n>0 


1 f e 2 b(n) 2 


2 \ n 


n! V 21 V 


+ E 


1 ( £ 2 b(ri) 2 \ "+ 1 / 2 


n> 0 


n! V 21 V 


from which we conclude that 
E(e' 


e\^{fn)-Vn{f n )h = ( 1 + fEE V - ( 

’ \ V2NJ^- ] 


n> 0 


2 \ n 


ra! V 2 N 


= 1 + 


eb(n)\ 

V2Nj 


^/2N)b{nf 


We end the proof of the lemma by replacing e by e\/~N. □ 
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